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Abstract 



\ We study the propagation of non-strange and strange meson modes in hadronic 
j> • matter considering both isospin and strangeness mixings induced by quantum fluc- 
. tuations in the medium. Baryons are described using the Quark Meson Couphng 
[ model extended to include interactions of strange quarks. In particular we evalu- 
ff^ ■ ate the dependence of the meson masses on the baryonic density, the strangeness 
fraction and the isospin asymmetry of the medium. We have found a considerable 
admixture of strangeness and isospin in the cr-mode in the high density regime. 
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The study of the meson properties in a hot/dense hadronic medium is at present an active field of research, since 
it is related to fundamental aspects of the low energy regime of the strong interactions such as the formation 
of a plasma of quarks and gluons and the restoration of the chiral symmetry. A prominent role is played by 
the lightest vector mesons, as they are closely connected to the emission rate of dileptons in heavy ion collision 
experiments. On the other hand the scalar mesons still keep valuable information about the quark structure of 
hadrons as can be seen, for example, in the dubious composition of the almost degenerate oo(980) — /o(980) 
mesons. 

There exist a profuse quantity of worthful studies of the in-medium modification of meson properties, some of 
them are based on pure quark models such as the Nambu Jona-Lasinio or use descriptions in terms of hadrons 
only [1,2,3,4,5,6,7]. There were also some efforts to reconcile both aspects by using hybrid models such as the 
Quark Meson Coupling (QMC) [8]. 

A very interesting feature appearing in these treatments is the mixing effect, which combines states of different 
isospin or Lorentz components [4,5,6]. This dynamical violation of the lagrangian invariances is particularly 
important in considering meson propagation, because it opens new in-medium decaying channels. 
The aim of this work is to study qualitatively the modification of the meson properties in a dense medium as 
described by the baryonic collective modes in the relativistic random phase approximation (RRPA) . In order to 
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adjust to the conditions found in the experience, as for instance in relativistic heavy ion cohisions, wc consider 
matter with diverse strangeness fraction as well as isospin asymmetry. We take into account the quark structure 
of baryons by using the QMC model [8,9,10]. In our treatment we regard all the mixing polarizations allowed 
by the hadronic lagrangian proposed. 

The evaluation of the RRPA with both isospin and strangeness mixing effects is still lacking at present, and 
the qualitative description of the lightest hadron masses in such a wide variety of situations as found in the 
experience is also deficient. 

The QMC model relates the in-medium properties of hadrons with their quark structure. Baryons are rep- 
resented as non-overlapping spherical bags containing three valence quarks, while the bag radius changes 
dynamically with the medium density. The exchange of mesons coupled directly to the confined quarks provide 
the baryon-baryon interaction. 

To fit the model to our purposes we consider the strange /o(980) and ^ mesons as in Ref. [10], and the scalar 
iso- vector ao(980) together with the more commonly used a, uj, and p mesons. 

The mean field approximation (MFA) equation for a quark of flavor q (q = u, d, s), of current mass niq and /g 
third isospin component is 

{i^d, - m,* - gl^'^u - 5^7°^ - aVl tV)*^ = 0, (1) 



where we have used 



m„ = m, 



,-9la-glC-gpl 6. (2) 



In these Eqs. 6 and ( stand respectively for the ao and /o fields, the upper bar indicates mean values. The 

solution of Eq. (1) for a quark confined within a spherical bag of radius Rh, representing a baryon of class b, is 
well known and can be consulted for example in [8,9]. The baryon effective mass is given by 



where is the number of quarks of flavor q inside the bag and Qqb = [x^^ + (-Rfomg*)^]^/^. The eigenvalue Xqi, 
comes from the so called linear boundary condition at the bag surface [9]. The bag constant B^/^ = 210.86Mey 
stands for non-perturbative vacuum contributions, it is adjusted to get a proton bag radius Rp = 0.6 fm in 
vacuum. The zero-point motion parameters zqi, are fixed to reproduce the baryon spectrum at zero density. In 
our calculations we have taken m„ = = 5MeV and = IbOMeV for the current quark masses. 
The equilibrium condition for the bag radius is [11] 
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dkk^ 



(4) 



The sum on the r.h.s runs over all the baryonic species considered. 

We neglect the coupling between strange and non-strange quarks and mesons. To determine the remaining 
couplings we assume pure vector dominance together with S\J (3) flavor symmetry, which allows to regain the 



2 



SU (6) relations of the quark model, i.e. 

5^,0) = ^ns9a,uj^ 9s,p = l-^sl 9s,p^ 9c,4> = ^s9c,<p^ (5) 

where N^^ (N^) is the non-strange (strange) quark number inside the baryon b, with third isospin component 
I^, and in addition we have g'^ = —V^g^. Furthermore as the SU{2) flavor in the non-strange sector is almost 
perfectly realized, we take as identical the coupling of quarks u and d with every meson. 

Therefore we have five independent couplings g^, g^, gg, g^, and 5^, we fix their numerical values by fitting 
symmetric nuclear matter properties at saturation, namely baryonic density no = 0.15/m~^, binding energy 
Eb = —IQMeV and symmetry energy = 35MeV. The last condition gives a functional relationship between 
gg and g^, we have chosen an arbitrary point on this curve. To fix g^ we follow references [3,10,14] assuming 
a potential depth Us = AOMeV for the S hyperon embedded in symmetrical S matter at baryonic density 
ns = no. Thus we obtain ga''^ = 5.99, g^''^ = 3.00, gg''^ = 1.83, gp''^ = 4.50, and g^ = 4.48. 

The self-consistent equations for the meson mean fields are of the form 

X = -^E^x< (6) 

with = dMh* /dx, n'^ = if the generic meson x is a Lorentz scalar field, or = —5^, — S'x) ~9x^3' '^x ~ 
for the vector case. The scalar and baryon densities 



kp 



n" = —^Mb / dk — , =, n — 

TT^ J 







are functions of the Fermi momentum kpb- The meson quantum fluctuations will be considered later. For the 
hadronic masses in vacuum we have used the values mo- = 550, = 980, m^^ = 783, m,^ = 1020, ms = 984, 
nip = 770, Mp = Mn = 938.92, Ma = 1115.63, and Ms = 1318.11, expressed in MeV. 

The QMC model describes effective baryons propagating in a homogeneous background of classical meson fields 
as described above. Beyond this picture we assume a linear coupling between the mesonic fluctuations and the 
baryons emerging from the MFA, with coupling constants given in Eq. (5). The full meson propagator can be 
approximated in a non-perturbative approach by summing the ring diagrams to all orders, i.e. the RRPA. This 
is a common procedure in nuclear physics, and can be used to evaluate the density dependence of the meson 
masses among other properties. The QMC has been used previously to evaluate the vector meson masses, 
regarding bosons as bags on the same foot as the baryonic ones, but neglecting the quark structure of the a 
meson [8]. 

The one-loop proper polarization insertion describing the propagation between meson states a and /3 is given 

by 

na/j(g) = -iY.9i9'pj [Gb{q) Gbiq + k) F^] , (7) 
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with indexes a, (3 running over the whole set of mesons and its internal degrees of freedom. Further, Gb{p) 
stands for the baryon b propagator in the MFA, and the structure of Fq depends on the baryon-meson vertex, 
namely Tcr^ = 1; = r, F^,^^ = 7^, Fp = 7^ r. Eq. (7) contains ultraviolet divergences coming from the vacuum 
contribution which require an appropriate regularization. We adopt the dimensional regular ization procedure, 
preserving that Hapiq) = H^aig), and we shall discuss later the choice of the regularization points. Explicit 
expressions for Eq. (7) are similar to that given in [2,6], for instance. 

The formalism is best described within a generalized meson propagator in a matrix representation of dimension 
equal to the sum of the mesonic degrees of freedom. For example the free (7^°), and full {V) generalized meson 
propagators have respectively in its diagonal blocks the free and dressed meson propagators of the a, (, S, uj, (j), 
and p fields, and in the complementary spaces they have zeros or finite mixing polarizations arising from Eq. 
(7), respectively. The corresponding Dyson-Schwinger equation in matrix form can be used to solve for V{q): 

Piq)=[l-P\q)Uiq)]-'p\q) (8) 

The dielectric function e{q) = det [X — V^{q) n(g)] is defined so that its roots coincide with the poles of V{q). 
We study the temporal regime q = 0, where scalar and vector contributions to the equation e(q) = become 
independent; in addition the longitudinal and transversal branches of the vector mesons become degenerate. 
Hence we have two separate equations of the form: 

{D^' - Uxx)iD^' - ^f.t.)iD-' - n,,) - (d;^ - Haa) n^v -{D-^- n..) n^^ = o, (9) 

where the indices A, /Lt, v take definite meson labels, namely \ = 5, ji = a, and v = C, ioY the scalar case, and 
X = p, fx = Lo, and v = (f) for the vector case. 

Furthermore = q^ ~ Tria inverse free propagator for the a-meson. It must be noted that there is no 

S-(^ neither p-(p mixing, and this causes the splitting of the isovector mode in vacuum as well as in the isospin 
symmetric dense medium. This splitting appears because the a-5 and uj-p polarizations become null under 
these conditions, in view of the degeneracy assumed for the baryonic iso-multiplets. 
In order to extract finite values from the polarization we require at zero baryonic density [12] 

^ap{q' = Rip) = 0, 

^ a/3 



The last equation holds only for both a and /? corresponding to scalar mesons, where g'^ = for the isoscalars 
a = a,(, or g^ = g^Is for the isovector a = 6. The regularizing parameter is fixed at Lap = 10 [7,12]. Here 
Rap = Rpa are a set of regularization points that we choose in order to reproduce the physical meson masses 
at zero density. This requirement determines unambiguously Rss = ms, and Rpp = rup, otherwise we obtain 
a range of points (iJo-o-, ^o-C' ^Cc) (Rujuj: Ru^rf): Rtfxf))- The regularization points for the mixing polarizations 
including one iso-vector meson are not constrained, for the reasons mentioned in the preceding paragraph. In 
our calculations we have used R^a = n^-a, RaC, = RaS = 2.807, R(^(^ = 4.654, R^^^^ = 3.98, R^jfj) = R^^p = 6.09, 
R(j)(t> = 5.12 expressed in fm~^. We have numerically checked that the final conclusions are insensitive to the 
precise choice of these parameters. Instead, the choice of L has noticeable effects on the behavior of the a 
meson properties [7,12], and a detailed discussion within this context will be given elsewhere [13]. 
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To study the behavior of meson masses in a hadronic medium with strange and isospin content we have 
considered two situations: a) nucleons and A, b) nucleons, A and H in equilibrium against strong decay. In the 
last case the relative abundance of A and H is determined by the relation between their chemical potentials, a 
situation often treated in the literature [3,10,14]. These two different hadronic environments when considered 
at zero isospin asymmetry give practically identical results for the meson masses. Therefore we consider in the 
following only case (a). 

We take the total baryonic density n = (n" + + n^) running up to 4no, guarantying that baryonic bags do 
not overlap. For a fixed total baryonic density n we introduce the isospin asymmetry parameter t = (n" — nP) / n 
and the strangeness fraction S = n^/n, both S and t taking values in (0, 1). 

In Fig. 1 we plot the effective meson masses computed as the zeros of Eq. (9) , for t = and different strangeness 
fractions S. Although the meson quantum numbers are blurred in-medium due to the mixing effect, we can 
distinguish well defined masses branches continuously related to its zero density value. Therefore we keep the 

label of the meson which originates each branch. 

To explain the behavior of the meson masses shown in this figure one must take into account that the polariza- 
tion insertion is composed of two additive terms of opposite sign, namely the vacuum contribution II^^) and the 
density dependent part II^^^ which explicitly depends on the Fermi momenta of the baryons. Thus, neglecting 
mixing corrections one has ~ + n[.p -|- H^.^'', k = (^,5,uj,(j) and p. At low densities (n < no) vacuum 
fluctuations are dominant and tend to lessen the meson effective masses. However the density dependent part 
of the polarizations contributes in the opposite direction, attenuating or even reverting this trend at higher 
densities. The only exception is given by m* for which 11 (-^^^ raises at low densities due to the relatively small 
value of mo-. 

Within this approach the dependence on S for a fixed baryonic density can be explained by observing that 
n(^) is a negative decreasing function of S. On the other hand 11*^^^ decreases (increases) when S increases 
for the non-strange (strange) mesons, due to the depletion (population) of the baryons to which they couple. 

Furthermore the variation of 11*^^^ is stronger for the vector than for the scalar mesons. For instance, as S 
grows the total polarization of the branch increases due to the predominance of H^'^^ over H^^^ , therefore 

is enhanced. On the other hand for the branch, the growing of II^^^ is much weaker and II^^^ determines the 
decreasing behavior of with S. The numerical results shown below provide the justification for disregarding 
in first approach the mixing effects for the strange mesons. This reasoning docs not hold for the a branch 
because it carries a strong mixing, leading to a more involved behavior, but the final result is a decreasing 
dependence with S. 

On the other side, as can be appreciated from this figure, the lj — p mass difference for baryon densities around 
2no is close to the pion mass, increasing the probability of the p oj + tt decay. This fact in turn contributes to 
enhance the p meson width, in agreement with several theoretical predictions as well as with phenomenological 
observations. 

We can compare Fig. 1 with the results shown in Ref. [3] (Figs. 3 and 4 for g = IMeV) for the a,uj and ( 
mesons. It can be appreciated that the overall behavior with S agrees with our calculations within the consid- 
ered range of densities. 

The influence of the isospin asymmetry t over the meson modes at constant strangeness S = 0.5 is shown 
in Fig. 2. The qualitative description does not change substantially for other values of S. The strange meson 
masses are almost independent of t, as they are coupled only to the iso-singlet A. On the left hand panel, the 
(T-meson mass shows a slight enhancement with growing t, this effect turns to be appreciable for n > 2nQ. 
On the opposite panel of this figure, the isovector p-meson mass also exhibits a small enhancement with t for 
n > 1.5no. Interestingly the isoscalar a;-meson mass presents the inverse trend when t increases, therefore we 
conclude that for a fixed strangeness content, the splitting between the u> and p masses increases with the 
isospin asymmetry. These results can be compared, for example, with those of Ref. [4] where the same behavior 
of ml, p with t is found. 
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A measure of the mixing between two physical meson states a and /3 can be given by the mixing angle Q^p^ 
defined by [5] 



tan 2^a/3 = -TT^ I2 TT • (10) 



m 



Since our aim is to obtain a qualitative estimation of how much isospin and strangeness remain good quantum 
numbers for the collective modes, wc do not include the width of the original mesons in Eq. (10) in consistency 
with the RRPA. However quantitative calculations must include a imaginary contribution to the polarizations 
[4,15]. 

For t = the isoscalar-isovector mixing vanishes, remaining the isoscalar strange-non-strange mixing only. We 
have evaluated Eq.(lO) for a = a (cj), (3 = Q {(j)) at the pole pQ = m* (m* ) in terms of the baryonic density for 
various strangeness fractions S. The result is plotted on the upper left (upper right) panel of Fig. 3; as expected 
the mixing effect increases with n and S. We can observe that Oac^iqo = m*) becomes appreciable at n > no, 
and therefore the mixing effects in the a channel may not be neglected even at relatively low densities. In the 
case of 9ujij>{qo = ?7i^) it turns to be significative only for high densities (n > 3no). With respect to 9a(;i% = nT'^) 
and 6uj(j>{qo = rn*^) (not shown here), their amplitudes remain below 3.5" for all the densities studied. 
Because the strangeness mixing remains small for all but the a meson, we can estimate the isoscalar-isovector 
mixing for t 7^ by using Eq.(lO) also. In the lower left panel of Fig. 3 it can be seen that scalar a — S mesons 
mix only appreciably in the a branch at densities n > 3.5no. Instead, non-strange to — p vector mesons exhibit 
a noticeable isospin mixing of approximately equal magnitude in both branches, as shown in the lower right 
panel of this figure. However, inclusion of the u and p widths is expected to decrease these values [4,15]. By 
raising the strangeness fraction a sizeable decrease of the mixing angle is obtained for both scalar and vector 
mesons, in the considered range of n. 

We conclude that the in-medium propagation of the a-channel concentrates non-negligible admixtures of isospin 
as well as strangeness in the high density regime. The vector u and p modes are affected by isospin but not 
by strangeness mixing. On the other hand the S, ^, and (f) modes propagate in almost pure states for all the 
densities studied. 

The effects discussed for vector mesons could be observable for instance, in the low mass dilepton production 
from heavy ion collision. The in-medium modification of the vector mesons has been successfully invoked to 
explain the gross features of the dilepton spectra from the CERES collaboration [16], in the frame of the Vector 
Meson Dominance model. A key role is assigned to the in-medium dropping of the p-meson mass. In order to 
get this result at the tree level, a direct coupling of the a to the co/p fields was proposed [16]. In our approach we 
obtain the same qualitative behavior within the RPA, and also a description of the mixture of states. Because 
of the smallness of the uj — (p mixing shown in Fig. 3 and the fact that both mesons decay largely after the freeze 
out, we do not expect a significative modification in these channels. The uo — p mix is large enough at almost 
all densities to have significative consequences. Because the p and uj mesons decay mainly into two and three 
pions respectively, in the dense hadronic medium the alternative mechanism p — ^ a; — ^ Stt is viable, having the 
effect of removing part of the strength between the peaks corresponding to both mesons and transferring it to 
the left of the p peak in the dilepton spectra. 

Some points deserve to be investigated within this model, for example the meson coupling to Goldstone bosons, 
the fact that the non-overlapping bag hypothesis is violated in the extreme density regime used in astrophysical 
applications [11], and the medium dependence of the scalar coupling constants which play a significative role 
in chiral models [17,18]. 
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Fig. 1. The effective meson masses as a function of the density for t = and several strangeness fractions 
5 = 0, 0.25, 0.5, 0,75, 1. The left (right) panel corresponds to the scalar (vector) meson branches. The arrows 
indicate the direction of growing S. 
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Fig. 2. The effective meson masses as a function of the density for S = 0.5 and several isospin compositions 
t = 0, 0.5, 1. The left (right) panel corresponds to the scalar (vector) meson branches. The arrows indicate the 
direction of growing t. 
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Fig. 3. The mixing angles defined by Eq. (10) as a function of the density. The upper panels show strangeness 
mixing for f = and several strangeness fractions S = 0, 0.25, 0.5, 0, 75, 1. The upper left (upper right) panel 
corresponds to the scalar a-C (vector Lo-(f)) mix evaluated at the a (u) meson branch. The arrows indicate the 
direction of growing S. Lower panels show the isospin mixing angles for t = 0.5 and the strangeness fractions 
S = 0, 0.5. The lower left (lower right) panel corresponds to the scalar a-d (vector Lo-p) mix evaluated at both 
the a and 6 {lo and p) meson branches. 
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